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We investigate electronic conductivity in layered metals in magnetic field in the weakly incoherent 
limit, when the interlayer transfer integral is smaller than the Landau level broadening due to the 
impurity potential, but the interlayer electron tunnelling conserves the intralayer momentum, ft 
is shown that the impurity potential has much stronger effect in this regime, than in the quasi-2D 
metals in the coherent limit. The weakly incoherent regime has several new qualitative features, 
not found in the previous theoretical approaches. The background interlayer magnetoresistance in 
this regime monotonically grows with increasing of magnetic field perpendicular to the conducting 
layers. The effective electron mean free time is considerably shorter than in the coherent regime 
and decreases with magnetic field. This enhances the role of higher harmonics in the angular 
magnetoresistance oscillations and increases the Dingle temperature, which damps the magnetic 
quantum oscillations. 
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The crossover between coherent and incoherent electron transport in the layered metals attracts great attention, both 
theoretical and experimental. [§-IT^| for its influence on the properties of high-temperature cuprate superconductors, 

"£^)i|ganic metals, heterostructures, and many other layered materials. This crossover can be driven by temperature T, volume 
rimpurity concentration external magnetic field B = (B Xl B y , B z ). In magnetic field this crossover in conductivity is very 
'"["pronounced because it qualitatively changes the magnetoresistance behavior. 

""^J i The electronic conductivity in magnetic field is widely used to investigate the electronic structure of various metals. In 
Strongly anisotropic quasi-2D layered metals, when the interlayer transfer integral t z is much smaller than the Fermi energy 
G?p the influence of magnetic field on conductivity has many specific features. One has to separate several different regimes 

^interlayer magnetotransport, depending on the ratios of three energy parameters: the interlayer transfer integral t z , 
the inverse mean free time Tq — fi./2-ro due to the impurity scattering, and the Landau level (LL) separation hio c , where 

^~~d; c = eB z /m*c is the cyclotron frequency. 
^ When the interlayer transfer integral is larger than the Landau level separation, t z > Hlj c , the 3D electronic dispersion is 

^^vell defined and given in the tight-binding approximation by 

e 3D (k) ^e(k x ,k y )-2t z cos(k z d), (1) 

^^vhere e (k x , k y ) is the in- plane electron dispersion and d is the interlayer spacing. Then the classical magnetoresistance shows 

3~^amaji oscillations. flol fla| which are used to determine the in- plane Fermi momentum. The magnetic quantum oscillations 
MQO) of interlayer conductivity in this case have beats of amplitude. [l7j and these beats are shifted with respect to the 
. b^ats of MQO of magnetization or of the other thermodynamic quantities. [18l420l | The slow oscillations also appear in the 
interlayer conductivity, which can be used to separate relaxation times from different scattering mechanisms. [20l l2l| 
Kj* When the interlayer transfer integral is smaller than the Landau level separation, t z < hui c , the beats of MQO disappear. 
%_This limit happens in strong fields in very anisotropic metals. If the interlayer transfer integral is still larger than the LL 
Ciroadening, t z > T, the dispersion ([1]) survives, and the MQO can be described by the "coherent" theory in Refs. [Tcl - t^ . 
Note, that the LL broadening T is larger than To in strongly anisotropic metals close to a stack of isolated conducting layers 
[see Eq. (|20]) below]. 

In the very anisotropic dirty limit, when the interlayer transfer integral is the smallest parameter, t z < hu c ,T, the tradi- 
tional 3D approach fails to describe the interlayer magnetoresistance. For example, in this limit, the experimentally observed 
interlayer magnetoresistance grows with increasing of the out-of-plane magnetic field strength B not only in the maxima, 
but also in the minima of MQO (see, e.g., Refs. fsj-fiol fl2|). The angular dependence of the background magnetoresistance 
also has many unusual features in this regime. [§. 1 12| This change of the magnetoresistance behavior as the magnetic field 
strength or the impurity concentration increase was called the " coherence-to-incoherence crossover" . It has been observed 
in various compounds and attracted the considerable theoretical attention. [l|, H, Q The term "weakly incoherent" has been 
introduced[l[ to separate this regime from the coherent 3D limit t z > h/r, and from the completely incoherent regime, where 
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the electron tunnelling to the adjacent layers does not conserve the in-plane electron momentum. The completely incoherent 
interlayer electron tunnelling happens when it goes via resonance impurities, 0, @, H| due to interaction with phonons[3] 5] 
and in some other models. 

The theory of weakly incoherent magnetoresistance in Ref. [l| is based on the phenomenological Green function (see Eq. 
(53) of Ref. which is equivalent to 

n k ^ ^ n ^ 

Here j is the number of conducting layer, related to the z-coordinate as z = jd. The Landau level number n and the 
momentum component k y form the complete set of quantum numbers of the 2D electrons in magnetic field with free electron 
dispersion 

e„ = hu c [n + 1/2) . (3) 

In magnetic field B = (B x , 0,B Z ) (we may choose B y = without loss of generality because the in-plane dispersion is uniform) 
the electromagnetic potential in the Landau gauge is A — {zB y , xB z — zB x , 0). Then the 2D electron wave functions are 

* ntky j (x,y) =* n (x- l 2 Hz [k v +jd/l 2 Hx ]) e ik »«, (4) 



^ _ exp (^/2P Hz ) H n (x/l Hz ) 



where 



H n {x/Ihz) is the Hermite polynomial and, for brevity, we introduced the notation of magnetic length components Ih x = 
\Jhcj eB x and Ihz = eB z . In the Green function in Eq. ([2]) the impurity scattering produces only the imaginary part 

of the self-energy iTq, which is independent of the quantum numbers {n, k y , j} and of the magnetic field strength B. This 
approximation is incorrect in the weakly incoherent regime, as will be shown in Sec. II below. 

The Green function in Eq. ([2]) is not suitable to study the MQO, because the MQO of the electron density of state (DoS) 
at the Fermi level lead to the similar oscillations of the electron self-energy, which must be taken into account in the theory 
of MQO.[19L l20l. l22M24j For electrons with 3D dispersion, as in Eq. (p}, in the Born approximation and after averaging over 
the impurity configurations, the imaginary part of the self-energy is proportional to the density of states, i.e. it acquires the 
oscillating energy dependence: 

T = T(e)=T [l + p(e,B)/p ], (6) 
where p (e) and po are the electron DoS with and without magnetic field. The electron Green function 

G o (r r ie)= V ^, J ^^X fc , J ( a! i.yi)« i,C ' ( ' 1 -" ) (7) 
M 2 ' J ' J e-e 2D (n,k y ) + 2t z cos(k z d)-iT(e) ■ { ' > 

n,k y ,k z y 

with r (e) from Eq. ^ has been substituted to the Kubo formula in the calculation of MQO of interlayer conductivity a zz in 
quasi-2D metals in Refs. (l9l 12(1 [22| - [24l | The completely incoherent hopping mechanism of the interlayer magnetotransport, 
which does not conserve the in-plane electron momentum during the interlayer hopping, has also been suggested @ to explain 
the exponential growth in interlayer magnetoresistance with decreasing temperature. However, all these approaches cannot 
explain the monotonic increase of magnetoresistance with increasing B in the minima of MQO, observed in /3"-(BEDT- 
TTF) 2 SF 5 CH 2 CF 2 S0 3 .li,[3 

Below we reexamine the approach based on Eqs. ([2])-([7]), in the weakly incoherent limit hw c > Tq > t z . We argue that 
Eq. ^ does not hold in this limit, and derive the different formula for the Green function. Then we calculate the interlayer 
conductivity with the new Green function and show that the new result considerably differs from that in the "coherent" 
theory in Refs. fll. Ii~8l - [2~ij . This explains several new qualitative features of MQO and of the angular dependence of interlayer 
magnetoresistance observed in the weakly incoherent limit. 



II. THE MODEL 



The electron Hamiltonian in layered compounds with small interlayer coupling consists of the 3 terms 

H = H Q + H t + Hi. 



(8) 
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The first term Hq is the 2D free electron Hamiltonian summed over all layers: 



H = ^£2D {m)cZ ,<■. 



+ 



where {m} is the set of quantum numbers of electrons in magnetic field on a 2D conducting layer, Eid (fri) is the corresponding 
free electron dispersion given by Eq. ([3]), and c^(c m ) are the electron creation (annihilation) operators in the state {m}. 
The second term in Eq. @ gives the coherent electron tunnelling between two adjacent layers: 

H t = 2t z [ dxdy[*Ux,y)y j ^ 1 (x,y) + $]_ 1 (x,y)* j (x,y)}, (9) 



where ^j(x,y) and ^(x, y) are the creation (annihilation) operators of an electron on the layer j at the point (x,y). We 
call this interlayer tunnelling Hamiltonian " coherent" because it conserves the in- layer coordinate dependence of the electron 
wave function (in other words, it conserves the in-plane electron momentum) after the interlayer tunnelling. The last term 

Hi = Y,Vi{r) (10) 

i 

is the impurity potential. The impurities are taken to be point-like and randomly distributed on the layers. The impurity 
distributions on any two adjacent layers are uncorrelated. The potential Vi (r) of any impurity located at point is given 
by 

V t (r) = US 3 (r - r<) . (11) 

In the 3D limit, when the interlayer transfer integral t z is much larger than the electron level broadening T due to the 
impurity scattering, the impurity potential Hj can be considered as the small perturbation for the electrons with dispersion 
([T]). In the Born approximation this gives T — irriiU 2 p (Ep) in agreement with Eq. ([5]), where rij is the volume impurity 
concentration, and p (Ep) is the DoS at the Fermi level. This leads to the standard theory of magnetic quantum oscillations 
in Q2D metals. 17] [20] In the opposite limit, t z -C F, tiLo c , the interlayer hopping t z must be considered as a perturbation for 



the disordered uncoupled stack of 2D metallic layers, where Eq. ([6]) is no more valid. |25l - l31| 

The 2D metallic electron system in magnetic field in the point-like impurity potential has been extensively studied. [25U3"i| 
The point-like impurity potential leads to the broadening of the Landau levels, which is described by the density of states 
(DoS) distribution function D(E). Since each Landau level without disorder is strongly degenerate, even weak impurity 
potential lifts this degeneracy and leads to the considerable broadening of the Landau levels. The electron Green functions 
acquire a cut instead of the pole as in Eq. ([2]). In the self-consistent one-site approximation, the Green function is given by 

G ( ri > r2 ' £ ) = E *°X (»*X*, (ri)G (e, n) , (12) 

where 



G (E, n) = ^MljWBB, (13 ) 
V ' ' 2EE g ' V ' 

and the DoS D (E) = (— 1/tt) IuiGe (E) on each LL is described by the dome-like function [25jj 

, x \/{H — Ei) (E 2 — E) , s 

D ^ = ~ 2«\E\E S ■ (W) 

where the electron energy E is counted from the last occupied LL, E = e — S2D (n, k y ), and 

E g = V /2nl 2 Hz . 

Here Vq = U \ip (zi)\ 2 w U/d, is the 2D analogue of the strength U of the point-like impurity potential: 

Vi(x > y)=V 5(x-Xi)5(y-y i ) > (15) 
and ip (z) is the out-of-plane electron wave function. The boundaries of the DoS dome in Eq. (fl4|) are 



E x = E g (y/cl - l) 2 , E 2 = E g {^-+l)\ 



(16) 
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where Cj is the ratio of the 2D impurity concentration, N{ = riid, to the 2D DoS on one LL, Nll — l^vr/f^: 

Cl = N t /N LL = 2<irl 2 Hz nid. (17) 

The function in Eq. (fT4"f is normalized to unity, J D (E) dE = 1. D (E) converges at the point E = 0, because this point 
lies outside the DoS dome E\ < E < E%. The LL broadening 



L B = (E 2 -E l )/2 = 2E g JT l o, VB. (18) 

The ratio 



L B _ 2V ^/n l d/2nl 2 Hz _ 2U^ ni N LL /d _ 2 y/mlPp (E F ) hiu c _ [mT c 



L nniU 2 p(E F ) 7rniU 2 p(E F ) irn,U 2 p(E F ) V ttL 

grows as \f~B~ in high magnetic field. Eqs. (|18p . (TTT)l) give the correct asymptotic for the LL broadening in strong magnetic 
field. In weak magnetic field, when huj c <C Tq, the mean scattering time tb related to level broadening as Lb = S/2tb and 
entering the Drude formula, does not depend on the value of magnetic field along the conductivity: tb = Tn.[32j To get the 
correct asymptotic behavior for Lb both in strong magnetic field and at B — 0, one can take the simple function 



L 



(4^ c /ttLo) 2 + 1 



1/4 



(20) 



More realistic models of the finite-range impurity potential, and more accurate calculation of the DoS, including the many- 
site corrections, lead only to the small tails of the DoS dome. 26] [28] 31] The number of electron states in these tails is much 
less than the number of states in the DoS dome and can be neglected. However, to include these tails into account and to 
simplify the subsequent calculation, one can take the Lorentzian DoS distribution with the same broadening: 

D {E) ~ Tr^ + rl) n — • (21) 

Combining Eqs. (TT2"]). ([Ml) and we obtain 

G(n, T2, £ = V ' " 'S , (22) 
e - £„ - iT B 

This Green function will be used in the next section to calculate the interlayer conductivity. 



III. CALCULATION OF CONDUCTIVITY. 



The interlayer conductivity a zz , associated with the Hamiltonian ([9]), can be calculated using the Kubo formula and 
the formalism, developed for the metal-insulator-metal junctions[33j The conductivity is expressed via the electron on the 
adjacent layers: 

azz = W^T y {l Srd2r ' \ ^A(r,r\j,s)A(r',r,j + l,e)[-n' F (s)}y (23) 
where the electron spectral functions 

A(r,r',j,s) = i [G A (r,r',j,s) - G R (r,r',j,e)] , (24) 

and the advanced (retarded) Green's functions GA(R)(r,r',ji e) include interaction with impurities. The product of two 
spectral functions in Eq. (|2"3")l rewrites as 

II = A(r, r', j, s)A(r', r,j + 1, e) = G A (r, r',j, e)G R (r', r, j + 1, e)+ (25) 
+G R (r, r',j, e)G A (r', r, j + 1, e) - 2ReG R (r, r', j, e)G R (r', r,j + 1, e), 

In addition to the terms with the product of Green functions GaGr, the expression for conductivity also contains the term 
— KeG R G R , which becomes important when MQO are considered. |20l [22^ Eq. (46) and subsequent formulas in Ref. [l[, 
where this term is omitted, can be applied only when MQO are disregarded. In strong magnetic field especially in the 
layered metals, on contrary, the MQO are very strong. 
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The angular brackets in Eq. (|23[) mean averaging over impurity configurations. Since the impurity distributions on each 
layer is uncorrelated with other layers, one can perform this averaging separately for each spectral function independently, 
which gives: 



HL x L y 



</-/•</-,-' / ^{A(r,r',j,e)) {A(r',r,j + l,e))[-n' F (e)} 



(26) 



The averaged Green (or spectral) functions are translational invariant: (G R (r , r' , j , e)) = (Gn(f — r',j,E)). 

If the magnetic field is tilted by angle 9 with respect to the normal to the conducting planes, B = (B x ,0,B z 
(Bsm9,0,Bcos9), the Green functions on two adjacent layers acquire the phase shift [see Eq. (49) of Ref. [l[]: 



G R (r, r', j + !,£) = G R (r, r', j, e) exp {ie [A (r) - A (*•')] /h} , 



where 



Substituting Eq. (J27J) into Eq. ([25]) we obtain 



and 



A(r) = -yB x d = -yBdsinf 



n = 2G A {r, r',j, e)G R (r', r, j, e) cos {e [A (r) - A (r')] /ft} - 
-2Re [G R (r, r',j, e)G R {r', r, j, e) exp {-ie [A (r) - A (r')] /h}] 



2e 2 t 2 d 



Jd 2 r J ^[-n' F (e)]x 



x |KG K (r,e)>rco8 
- Re 



2 / e5z/d 



V ft 



■ sin ( 



/r , / u 2 / ieByd . 
(G R (r,e)) exp I — : — smfl 



The term in the third line of Eq. d^HJ) is absent in Eq. (50) of Ref. 
In the magnetic field perpendicular to the conducting layers 



2e 2 t 2 d 



\(G R (r,e))\ 2 -Re(G R (r,e)) 2 [-n' F {e)]. 



The integration over r for the Green function of the form (|12[) is very simple and gives 

2e 2 t 2 z dN LL f de 



With the approximate Green function, given by Eq. (f2"2")l . Eq. (|5T|) becomes 



2eHidN LL f de 



n 



2n 



E 



-n' F {e)]2T 2 B 



[{e-e n f + n 



(27) 



(28) 



(29) 



(30) 



(31) 



(32) 



The sum and integral in Eq. (|32l) is calculated in a standard way, transforming the sum over LL into the harmonic sum by 
applying the Poisson summation formula: |34j 



oo CO „oo 

E /(»)= E / ^ lkn f(n)dn 



(33) 



where a € (no — 1; no). Then, performing the integrations, we obtain 



<j zz =<j {B) ^2 (-l) fc exp 



k= — oo 



2tt (ikfj, - \k\T B ) 



ftiVr. 



2kT7 2 T/hu c 
sinh {2k-K 2 T/huj c 



2^|fc|r B 



(34) 



G 



where 



*o(B) = -^, 35 
vp = Nll/Hwc is the DoS at the Fermi level in the absence of magnetic field. Eq. would coincide with Eqs. (17)-(21) 



of Ref. [22| if To and T e in these equations are replaced by Tb- Note, that the nonoscillating part of conductivity <7o (B) 
is now a function of magnetic field, because oc \J~B in strong magnetic field. This observation contradicts the previous 
theoretical results fl| (22j [23j . also developed for the almost 2D case, because in these papers the LL width Tb was incorrectly 
taken to have no monotonic dependence on magnetic field, Tb = Tq + T, where T rapidly oscillates around zero. 

Let us now compare how strongly our result differs from the previous results [see, e.g. Eqs. (17-21) of Ref. [22|]. In Fig. 
[T]we plot the MQO of resistivity R zz (B) = l/a zz , calculated using Eq. (jMf with Tb given by Eq. (|2T)]) [solid blue line] and 
with Tb = Tq [dashed red line]. The difference is evident: the interlayer magnetoresistance shows monotonic growth and 
weaker oscillating amplitude with increasing magnetic field, than in the old result. 




FIG. 1: The MQO of resistivity R zz (B) — l/a zz , calculated using Eq. with Fb given by Eq. pp) [solid blue line] and with 

r_B = F) [dashed red line]. 

In tilted magnetic field the calculation of Eq. (j2"5)) performed in Ref. ]l| can be applied with the new magnetic-field- 
dependent value T B instead of Tq, which gives [compare to Eq. (1) of Ref. [1|] 

= (B) I [Jo (k)} 2 + 2 J2 - |_ J ; ( " )]2 ) , (36) 

where k = kpdt&nO and 

T B = h/2T B = r (r /r B ) . (37) 

There are two differences between this formula and Eq. (1) of Ref. [l|. First, the higher harmonics in AMRO are weaker- 
damped in Eq. ([31))) because of smaller value of tb oc 1/vB, Second, as we noted before, the background conductivity 
<t (B), given by Eq. (|3"5"1) . decreases as \j\fB in strong field. 

The higher harmonic in Eq. (|3!)|) play considerable role in AMRO. To illustrate this, in Figs. 12131 we plot the angular 
dependence of conductivity a zz (9) given by Eq. (|5rJ|) with t b — h/2T B = tq {Tq/T b ) and t b = To. For simplicity, we take 
the axially symmetric case, i.e. the symmetric in plane electron dispersion. One can see from Figs. 12131 that in the minima 
of conductivity, i.e. at the Yamaji angles, the replacement To — > t b is very important. The predicted value of conductivity 
at the Yamaji angles with tb given by Eq. (|37|l is much larger than with tb — tq. This difference increases with increasing 
of magnetic field. The positions of the conductivity minima, i.e. the Yamaji angles, also slightly shift after the replacement 
t T B in Eq. {36|) [see Figs. EE]. For the first Yamaji angle at B = 5T (see Fig. [2]) this shift M Yam w 1.7°. 



IV. DISCUSSION 



Let us formulate the main difference of the present approach to the calculation of interlayer conductivity in the weakly 
incoherent regime compared to the previous methods, developed in Refs. [19] 22] [20] [23] 24| to calculate the MQO of conduc- 
tivity. In these papers the impurity potential is considered as a small perturbation on the background of a free electron gas 
with well-defined 3D electron dispersion given by Eq. (TI]). Hence, the impurity scattering was taken into account only by the 
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FIG. 2: The angular dependence of conductivity o zz , calculated using Eq. (|36|l with tb given by Eq. Q37|) [solid blue line] and with 
tb = to [dotted red line]. The dashed green line gives the difference between these two curves. The parameters for this plot are 
k F d = 2, m* = 2m e , B = 5T, T = IK. 
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FIG. 3: The angular dependence of conductivity a zz , calculated using Eq. (|36|) with tb given by Eq. (|37[) [solid blue line] and with 
tb ~ tq [dotted red line]. The dashed green line gives the difference between these two curves. The parameters for this plot are the 
same as in Fig. [2] besides the larger value of magnetic field B — 20T. One can see that in the minima of conductivity, i.e. at the 
Yamaji angles, the difference is very strong. The predicted value of conductivity at the Yamaji angles with tb given by Eq. (|37[) is 
much larger than with tb = tq. 

imaginary part of the electron self-energy given by Eq. (JBJ) . Even less accurately the impurities are treated in Ref. [l| , where 
the constant electron mean-free time has been used to include the interaction with impurities. This is correct only in the 
coherent limit, when the interlayer transfer integral is much larger than the LL broadening, and the electrons, moving in a 3D 
metal, are scattered by impurities. In the weakly incoherent regime, when t z < T, this is incorrect, because for a 2D electron 
system in magnetic field the impurity potential has much stronger effect than in 3D. Simply, in a 3D electron system the 
electrons after scattering by an impurity move away in the interlayer direction and never return to this impurity. Therefore, 
this impurity only leads to the single scattering of this electron into some other state, which is well described by the constant 
electron mean- free time r , or equivalently, by the constant imaginary part Tq of the electron self-energy. In 2D electron 
system in magnetic field, the electrons after scattering return to the same impurity after the cyclotron period. Therefore, the 
impurity has permanent influence on the electron state, considerably shifting the electron energy and modifying the electron 
states. Hence, in the weakly incoherent regime, when t z < T, the interlayer hopping term ([9]) in the Hamiltonian jSJ, rather 
the impurity potential (1101) . must be considered as a small perturbation. Therefore, to calculate the interlayer conductivity, 
we start from the stack of isolated 2D disordered conducting layers in magnetic field, where the effect of impurity potential 
is considered much more accurately, at least in the self-consistent one-site approximation. Then we substitute the obtained 
electron Green functions to the Kubo formula for the tunnelling conductivity between adjacent conducting layers. The effect 
of impurities in the final results turned out to be much stronger than in the previous approaches. Phenomenologically, this 
difference can be taken into account by the replacement of the initial level broadening Tq by the new value given by Eq. (JBJ) . 

One can also obtain Eq. (|2"2l with the new value of T given by Eq. (|2U)) using different arguments. The physical origin 
of large DoS broadening in Eq. (|14[) is not the finite lifetime r of electron states, with is mathematically described by the 
imaginary part of the self-energy ImS = Tq = fo/2r, as in the 3D limit. On the 2D layers the LL broadening comes from the 
energy shift of each electron state, which is described by the state-dependent real part of the electron self-energy Re£. The 
averaging of the electron Green function in Eq. @ over the impurity configurations is independent on each conducting layer, 
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since the impurity distribution is assumed to be uncorrelated. Then, the coordinate part of the Green function remains of 
the form ([2]) with the bare electron wave functions in numerator [see Eq. (| 12[) and Appendix] , but the denominator acquires 
the real part of electron self energy, which is distributed with the DoS function D (E) : 



(G&(ri,r2,i,e)) = / dED(E) 



^ e-E-hw c {n + 1/2) - ir ' 



The triangular brackets indicate averaging over impurity configurations. Substituting the approximate Lorentzian DoS 
distribution, given by Eq. (I21|) . one can easily perform the integration over E and obtain 

(G R {r^ w) ) « g e -^ e ( n + i /2 )- i( r + r fl ) - (38) 

This Green function differs from Eq. by the increase of the imaginary self-energy part: To — > To + T b with Tb given by 
Eq. ([15]). This is almost equivalent to Eq. (J22J with T B given by Eq. (|20]). 

Unfortunately, the proposed analysis considers only the limiting case Tb ^> t z , when is given by Eqs. (TTS)) or (|2"0)l . but 
it is not accurate for the intermediate case Tb ~ t z , where the crossover from the coherent to the weakly incoherent regime 
takes place. The phenomenological formula (|20[) gives only a qualitative dependence Tb {B z ) in this region. The crossover 
from the coherent to the weakly incoherent regime may be driven by the disorder (impurity concentration) or by magnetic 
field B z . The latter happens, because with the increase of magnetic field the effective LL broadening Tb also increases 
[see Eq. ([T8"]) ] and at some crossover field B c ~ t^m*c/roefi. it becomes greater than the interlayer transfer integral t z . To 
calculate the exact value B c of the crossover field and to describe the behavior of interlayer conductivity in this region one 
needs to calculate the electron Green function in layered metals with impurities and magnetic field in the crossover region 
Tb ~ t z . This is an interesting and still open problem. 

Above, we have not studied the MQO in the tilted magnetic field. The second term in the curly brackets in Eq. (|2"5)) does 
not contribute to the background magnetoresistance, but it affects the MQO. This term amplifies the MQO and modifies 
the angular dependence of the MQO amplitude. This modification is a fine effect which is harder to measure. The angular 
dependence of MQO amplitude is also affected by the Zeeman splitting and possible magnetic ordering. 

If the normalized point-like impurity concentration a < 1, the Nll — Ni electron states on each LL left degenerate, and 



besides the DoS dome the sharp S (E) term in the DoS survives. [2 7j However, as has been shown in Ref. [30(, the numerous 
weak impurities and the impurities, situated far from the conducting layers, are important for the lifting of the LL degeneracy 
in all layered materials. For achievable magnetic field even in the pulsed magnets B < 100T, Ih z > 25A Therefore, the 
typical normalized impurity concentration is greater than unity, Cj > 1 , and one can use the one-maximum DoS distribution 
as in Eq. (TPI|) , 

Above we have shown that the weakly incoherent regime strongly differs from the coherent limit. It also differs from 
the completely incoherent limit, where the new mechanisms of the interlayer electron transport, including the electron 
interlayer transport via resonance impurities 0, 0, HI and the hopping conductivity between completely localized states @, 
play important role. One difference of the weakly incoherent regime from the completely incoherent one is that the angular 
magnetoresistance oscillations (AMRO) are not damped, being of the same amplitude as in the coherent regime. Only higher 
harmonics in AMRO increase, making AMRO maxima less pronounced. The second difference is that the temperature 
dependence of conductivity in the weakly incoherent regime is the same, as in the coherent limit (usually, metallic), while the 
temperature dependence of the hopping conductivity [6| is exponential. Therefore, the weakly incoherent regime of interlayer 
magnetotransport is a separate regime, which should be distinguished from the coherent and completely incoherent limits. 



V. SUMMARY 



We reexamine theoretically the conducting properties of layered metals in the "weakly incoherent" regime, when the 
interlayer transfer integral t z is much less than the Landau level broadening F^ due to the interaction with impurities. 
The magnetic quantum oscillations and the angular dependence of interlayer conductivity in this regime are calculated. We 
obtain that both these effects in the weakly incoherent limit considerably differ from the coherent regime. This contradicts 
the previous theoretical results. [l[ [22[ The background interlayer conductivity a zz decreases with the increase of magnetic 
field B z according to Eq. (|35[) with Tb approximately given by Eq. (|20p . while in the coherent limit it remains constant 
(see Fig. [T]for illustration). The Dingle temperature of MQO also increases with magnetic field cx Tb- Meantime, in the 
weakly incoherent regime the angular oscillations of background magnetoresistance are not damped as in the completely 
incoherent mechanisms of the interlayer electron transport, considered, e.g., in Ref. [2], 6], [8]. On contrary, the damping of 
higher harmonics in the angular magnetoresistance oscillations is weaker than in the coherent regime [see Eqs. (|36p . (|37[) ]. 
This leads to the different picture of AMRO (see Figs. I2|.'>[) . Phenomenologically, the differences between the coherent and 
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FIG. 4: The first diagram for the electron self-energy with the intersection of the impurity lines. 

a a a 

f A A 

I / \ / 1 \ 

S=j+^ \ + / I \ + 

FIG. 5: The set of diagrams for the irreducible self-energy, corresponding to the self-consistent one-site approximation. The double 
solid line symbolizes the exact electron Green's function. 

weakly incoherent regimes can be taken into account by the replacement of the electron mean free time tq by the new value 



t~b ~ t~o/ (8w c To/7r) 2 + 1 in all formulas for MQO and for the angular oscillations of interlayer conductivity. 
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Appendix A: The in-plane electron Green function in the impurity potential 

Consider the noninteracting 2D electron gas in the potential of randomly distributed point-like impurity, as given in Eq. 
(|15p . The peculiarity of the two-dimensional electron gas in strong magnetic field in the presence of impurities is that the 
Born approximation of the scattering amplitude on each impurity is insufficient to describe the system. Physically, this 
means that an electron scatters many times by one impurity, because in magnetic field the electrons periodically return 
to the same point after passing along the cyclotron orbit. In the 3D case the diagram in Fig. [4] with the intersections of 
impurity lines is small by the parameter m/n e , where rii and n e are the volume impurity and electron concentrations. In 2D 
case in magnetic field there is no general proof that the diagrams with intersections of impurity lines are small. However, 
the calculations of the DoS in Refs. [HI, HH, [3l[ show that these diagrams only lead to the small tails of the DoS. Therefore, 
in the our subsequent analysis we keep only the diagrams without intersections of impurity lines. 

Now we proof by the method of mathematical induction, that, if one neglects the diagrams with the intersection of impurity 
lines, the electron Green function, averaged over impurity configurations, has the form of Eq. (|12j) with 

G(e,n) = 1/ [e - e„ - E„ (e)] . (Al) 

The energy e n of the n-th LL is given by Eq. and the electron wave functions k (r) by Eq. The self energy part 
£„ (e) for the n-th LL must be determined self-consistently and is given by the set of diagrams, shown in Fig. [5] In the 
self-consistent one-site approximation the self-energy part is 



1/4 



_ E - E g (1 - a) y/(E - gO (E 2 - E) 

n( ) 2 + 2 

The restriction given by Eq. (j!2[) is nontrivial because G (e, n) does not depend on k y . 

Without impurities, i.e. in the zeroth order of mathematical induction, Eq. ([12")) holds by definition. Assume, it holds for 
an arbitrary number N of impurities in the electron Green function Gn(t\, r%,e). When we add one more impurity center, 
the Green function Gjv+i {t\ , r% , e) is given by 

G N+1 (r 1 ,r 2 ,e) = J dr a G N {r 1 ,r a ,e)G N (r a ,r 2 ,e)T, (e,r a ) , (A2) 
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where £ (e) is given by the set of diagrams in Fig. [5] with the double-line standing for Gjv(r Q , r a , e). Performing the 
integration over k y in Eq. (1121) . we find 

G N {r a ,r a ,s) = ~ — : % 7-T- 

^ E - E n - 2^N,n {£) 

Therefore, £ (e, r a ) — S (e), and substituting Eq. ([12")) into Eq. (|A2I) and integrating over r a , we obtain 

GW+i(ri,r 2 ,e) = 2^—^ + c, / rfr a ^ — » ^ — » S(e) + .. 

7T ? rC y TL ^K/y 7i ? y 

^ e - e« - Sat,™ (e) ^ v e - e„ - £;v >n (e)J ^ e - e„ - S^ +1 ,„ (e) ' 

where 

Sjv+i,n (e) = S w ,n (£) + c»S (e) ■ 
Eq. (|A3j) has the form fj X2|) . which proves our statement. 



[1] P. Moses and R.H. McKenzie, Phys. Rev. B 60, 7998 (1999). 
[2] A. A. Abrikosov, Physica C 317-318, 154 (1999). 

[3] Urban Lundin and Ross H. McKenzie, Phys. Rev. B 68, 081101(R) (2003). 
[4] T. Osada, K. Kabayashi and E. Ohmichi, Synth. Met. 135-136, 653 (2003). 
[5] A. F. Ho and A. J. Schofield, Phys. Rev. B 71, 045101 (2005). 
[6] V. M. Gvozdikov,Phys. Rev. B 76, 235125 (2007). 

[7] D. B. Gutman and D. L. Maslov, Phys. Rev. Lett. 99, 196602 (2007) ; Phys. Rev. B 77, 035115 (2008). 
[8] M. V. Kartsovnik, P. D. Grigoriev, W. Biberacher, and N. D. Kushch, Phys. Rev. B 79, 165120 (2009). 

[9] F. Zuo, X. Su, P. Zhang, J. S. Brooks, J. Wosnitza, J. A. Schlueter, Jack M. Williams, P. G. Nixon, R. W. Winter, and G. L. 

Card, Phys. Rev. B 60, 6296 (1999). 
[10] J. Wosnitza, J. Hagel, J. S. Quails, J. S. Brooks, E. Balthes, D. Schweitzer, J. A. Schlueter, U. Geiser, J. Mohtasham, R. W. 

Winter, et al, Phys. Rev. B 65, 180506(R) (2002). 
[11] T. Valla, P. D. Johnson, Z. Yusof, B.Wells, Q. Li, S. M. Loureiro, R. J. Cava, M. Mikami, Y. Mori, M. Yoshimura, and T. Sasaki, 

Nature London 417, 627 (2002). 

[12] M.V. Kartsovnik, D. Andres, S.V. Simonov, W. Biberacher, I. Sheikin, N. D. Kushch, and H. Miiller, Phys. Rev. Lett. 96, 166601 
(2006). 

[13] J. G. Analytis, A. Ardavan, S. J. Blundell, R. L. Owen, E. F. Garman, C. Jeynes, and B. J. Powell, Phys. Rev. Lett. 96,177002 
(2006). 

[14] V. N. Zverev, A. I. Manakov, S. S. Khasanov, R. P. Shibaeva, N. D. Kushch, A. V. Kazakova, L. I. Buravov, E. B. Yagubskii, and 

E. Canadell, Phys. Rev. B 74, 104504 (2006). 
[15] K. Yamaji, J. Phys. Soc. Jpn. 58, 1520 (1989). 

[16] R. Yagi, Y. Iye, T. Osada, S. Kagoshima, J. Phys. Soc. Jpn. 59, 3069 (1990). 

[17] Shoenberg D. "Magnetic oscillations in metals", Cambridge University Press 1984. 

[18] PD. Grigoriev, M.V. Kartsovnik, W. Biberacher, N.D. Kushch, P. Wy der, Phys . Rev. B 65, 60403(R) (2002). 
[19] PD. Grigoriev, M.V. Kartsovnik, W. Biberacher, P. Wyder, [cond-mat/0108352| (2001) [unpublished]. 
[20] PD. Grigoriev, Phys. Rev. B 67, 144401 (2003). 

[21] M.V. Kartsovnik, P.D. Grigoriev, W. Biberacher, N.D. Kushch, P. Wyder, Phys. Rev. Lett. 89, 126802 (2002). 

[22] T. Champel and V. P. Mineev, Phys. Rev. B 66, 195111 (2002). 

[23] V. M. Gvozdikov, Phys. Rev. B 70, 085113 (2004). 

[24] T. Champel and V. P. Mineev, Phys. Rev. B 74, 247101 (2006). 

[25] Tsunea Ando, J. Phys. Soc. Jpn. 36, 1521 (1974). 

[26] Tsunea Ando, J. Phys. Soc. Jpn. 37, 622 (1974). 

[27] E.M. Baskin, L.N. Magarill, and M.V. Entin, Sov. Phys. JETP 48, 365 (1978). 
[28] E. Brezin, D.I. Gross, C. Itzykson. Nucl. Phys. B 235, 24 (1984). 

[29] The Quantum Hall Effect, Ed. by R. Prange and S. M. Girvin (Springer, New York), 1987. 

[30] A.M. Dyugaev, P.D. Grigor'ev, Yu.N. Ovchinnikov, JETP Letters 78, 148 (2003). 

[31] I.S. Burmistrov, M.A. Skvortsov, JETP Lett. 78, 156 (2003). 

[32] A. A. Abrikosov, Fundamentals of the theory of metals, North-Holland, 1988. 

[33] G. Mahan "Many-Particle Physics", 2nd ed., Plenum Press, New York, 1990 

[34] "Standard Mathematical Tables and Formulae", CRC Press, (1996) (or another similar book). 



